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1. Introduction 

Discovered more than 20 years ago, in 1995, by the CDF and DO collaborations |[T], ^ at 
the Fermilab Tevatron proton antiproton collider, the top quark is the heaviest elementary particle 
known today. Its high mass might indicate that it plays a special role in electroweak symmetry 
breaking and makes it a promising candidate as window to new physics. The detailed study of the 
top quark and its production mechanism is therefore crucial. 

In this article, measurements of top-antitop quark production {tt) cross sections, performed by 
CDF and DO in Run II of the Tevatron, are discussed. In Run II, lasting from 2001 to 2011, the 
collision energy was 1.96 TeV. The measurements of inclusive tt cross sections and their compar¬ 
ison to predictions from quantum chromodynamics (QCD) calculations allows to test the standard 
model (SM) in detail. Many new physics models with final states similar to that of tt production 
would change the measured cross sections significantly. Furthermore, differential measurements 
of tt cross sections are reported, which serve as test for perturbative QCD (pQCD) predictions, can 
be used to tune Monte Carlo simulations, and are a useful tool to look for new physics that would 
change the shape of kinematic or topological variables relative to the SM prediction. 

2. Inclusive Cross Section Measurements 

The production of tt pairs occurs via the strong interaction. At the Tevatron, approximately 
85% of the tt production happen via quark-antiquark annihilation, and 15% via gluon-gluon fusion. 
These proportions are approximately inverted at the LHC. This means that the sensitivity to beyond 
the SM (BSM) physics scenarios is different at Tevatron and LHC. 

The measurement of the tt cross section is done separately in different final slates. In fhe 
SM, fhe fop quark decays almosf 100% of fhe lime info a W boson and a b quark, wifh fhe W 
eifher decaying info a quark-anliquark pair, or a charged lepfon and fhe corresponding neulrino. 
The even! selecfion in each decay channel is adjusted fo fhe respective final sfale objecfs, wifh 
fhe goal fo enhance fhe dala sample wifh tt evenls and reduce fhe background, while keeping a 
reasonable sfafisfics. The decay channels considered are usually fhe dilepfon final slate, where 
bolh W bosons decay info eleclrons, muons or leplonically decaying fans, fhe semileplonic final 
sfale, wifh one W boson decaying info eleclron, muon or leplonically decaying lau, and fhe olher W 
boson info quarks, and fhe fully hadronic final slate, in which bolh W bosons decay info quarks. In 
fhe dilepfon final slate a low background can be achieved, bul fhe branching fracfion is small, while 
allhadronic evenls represenf a large branching fraction, bul a large background. The besl mixture 
for most analyses is usually the semileptonic channel, which mixes a good branching fraction with 
a manageable background. Depending on the channel, different approaches to the measurement of 
the cross section are done. In very clean channels usually a counting-method is sufficient, while 
less clean ones require the use of multivariate analysis techniques, in which several variables with 
small signal-to-background separation are combined into one discriminant. Given the two fj-Jets in 
the tt final sfale, a useful fool offen applied is fj-lagging, where properties of displaced fracks and 
secondary verfices from fhe B hadron decay are considered fo idenlify jels from b fragmenfalion. 

In Ihis article, analyses in Ihese Ihree main final sfales are presented. The main background 
conlribulion in dileplonic evenls comes from Z-f-jels processes, which can be simulaled using Monle 
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Carlo (MC) generators. In semileptonic events, the handling of the largest background contribution, 
Ik+jets, is done using simulated events and normalizing the yield using a data-driven approach. For 
multijet events, the largest contribution comes from QCD multijet events, which require data-driven 
modelling of the background. Further, smaller background contributions are single top, diboson 
and fake events, where the latter are events coming from jets misidentified as leptons in the detector. 

2.1 Tevatron Combination 

In 2014, CDF and DO published the first combination of the tt cross section from both col¬ 
laborations, using data samples between 2.9 and 8.8 fb^^ In particular, four analyses from CDF 
and two from DO were used in the combination. CDF provided one analysis in the dilepton fi¬ 
nal sfafe, where fhe number of evenfs with at least one f^-tagged jet was counted, two analyses in 
semileptonic events, where the number of events with at least one Zj-tagged jet was counted or a 
Neural-Network discriminant was built based on kinematic variables. The fourth analysis, in 
the allhadronic final sfafe, used evenfs wifh exacfly one or af leasf one Zj-fagged jef, performing a 
maximum likelihood fif fo fhe reconsfrucfed fop quark mass. The DO collaboration provided an 
analysis in fhe dilepfonic final sfafe, where a likelihood fif fo a discriminanf based on a Neural- 
Nefwork Z^-fagging algorifhm, was performed, and an analysis in fhe semilepfonic final sfafe was 
performed, in which evenfs wifh fhree and more fhan fhree jefs were splif info evenfs wifh zero, 
one or af leasf fwo Z^-fagged jefs. For signal dominated sub-channels a simple counting mefhod was 
fhen used, while in background-dominafed sub-channels a random fores! discriminanf was applied. 

The CDF analyses were combined using BLUE [Q], while fhe DO analyses were filled simiil- 
laneously wifh a likelihood til, where fhe systematic uncerlainlies were Irealed as nuisance param¬ 
eters. These fwo experimenl-specific combinations were fhen combined in a further step using 
BLUE, resulting in a ft cross section of rt;? = 7.60±0.41(slaf-|-sysf) pb |§] for a fop quark mass of 
172.5 GeV. Thisresullis compatible wifh fhe SMprediction of Op = 7.351q 27(scale-|-pdf) pb, cal¬ 
culated af nexf-fo-nexf-fo-leading order (NNLO) precision in QCD, with soft-gluon resummation 
to next-to-next-to-leading logarithmic accuracy [^]. The cross section measurements in different 
channels are compatible with each other, showing no evidence for possible BSM contributions. 
Eigure [I] shows the tt cross sections in the different channels by the two experiments, as well as the 
combination. 

2.2 New DO result 

Using the full Tevatron data sample of 9.7 fb^^, the DO collaboration recently updated their 
inclusive measurement of the tt production cross section in the dileptonic and semileptonic final 
sfafes. The analysis sfrafegy is similar fo the one used for the Tevatron combination: in dileptonic 
events, the discriminant trained for Zr-jet identification via a multivariate analysis technique is used 
to discriminate tt signal from the background. A template fit is performed to extract the cross 
section value. In semileptonic events, a topological discriminant, based on boosted decision trees 
(BDTs), is trained separately for events with two, three or at least four jets, and the cross section 
is extracted by simultaneously performing a maximum likelihood fit on the discriminants in these 
three sub-channels. The resulting cross section is Otj = 7.73 ±0.13(stat) ±0.55(syst) pb [^, in 
good agreement with the SM prediction. The main contribution to the systematic uncertainties 
comes from modelling of the hadronization. 
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Figure 1: Tevatron combination of the tt cross section by the CDF and DO collaborations [||]. 


3. Differential Cross Section Measurements 

Besides performing measurements of inclusive cross sections, the determination of the tt cross 
section differentially as function of various variables is important. Differential tt cross section mea¬ 
surements provide a test of pQCD, help to tune the simulation of tt events, and can yield additional 
insight into potential hints for new physics. For example, deviations in the tt cross section as func¬ 
tion of the invariant tt mass, could yield insight into the existence of heavy resonances decaying 
into a pair of top quarks, or possible new physics contributions to the tt forward-backward asymme¬ 
try can be probed by measuring differential distributions in variables related to the pseudorapidity 
of the top quark. 

Using the full Run II data sample of 9.7 fb^\ DO recently performed a measurement of the 
tt cross section as function of three variables: nitj, the absolute value of the rapidity of the top 
quark, and the transverse momentum of the top quark, pj^ ||^. The analysis was performed 
in the semileptonic final state, with at least one jet required to be Zi-tagged. The calculation of 
the variables \ \ and pj^ requires the assignment of final sfafe objecfs fo originate from fhe fop 

or fhe anfifop quark. For fhis assignmenf, a consfrained kinematic reconsfrucfion algorifhm is 
used, in which experimenfal resolutions are faken info accounf. In fhis algorifhm, fhe fop quark 
mass and W boson mass are fixed fo fheir known values, allowing fo defermine fhe z-componenf 
of fhe neufrino momenfum from fhe W boson mass consfrainf and reducing fhe number of possible 
jef-quark assignmenfs via fhe fop quark mass consfrainf. The solufion wifh fhe besf is taken for 
the construction of the top quark vectors. 

The measurement is defined for parfon-level fop quarks, including off-shell effecfs. There¬ 
fore, a correcfion for defecfor and accepfance effecfs has fo be done. A regularized unfolding, 
implemented in TUNFOLD } } ~\citelunfold is used for fhis purpose, where fhe regularization is 
based on fhe derivative of fhe disfribufions. To keep as much informalion as possible, fhe number 
of bins for fhe reconsfrucfed disfribufions is kepf fwice as high as fhe number of bins on parfon 
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Figure 2: The tt cross section measured as function of \}/°^\ (top), and pj’’ (bottom). The left figures show 
the cross section, the right figures the ratio of the measured cross section with respect to the approximate 
NNLO calculation 0. 


level. The contribution of different sources of systematic uncertainties are evaluated by changing 
the migration matrix and the background contribution. The largest uncertainties arise for high bins 
in nitj, and pj^. In Figure^ the measured tt cross section differentially in and pj^ is 
shown, together with predictions from Monte Carlo simulation and predictions from approximate 
NNLO calculations. In general, the agreement is good relative to the generator and approximate 
NNLO predictions. For ALPGEN [|TT1], the absolute normalization is too low, while in the \y‘°P \ 
variable the description by the MC@NLO generator [1^] is better than that of the approximate 
NNLO prediction. 

The tt forward-backward asymmetry has been measured by both CDL and DO in recent years 
to be somewhat larger than predicted by the SM [13]. In many BSM models, for example t! 
or axigluon models, this asymmetry would be associated to a change in differential distributions. 
Ligure ^ shows the differential tt cross section as function of compared to different axigluon 
and a t! model. The measurement excludes several of these models. 


4. Conclusion and Outlook 

The measurements of the inclusive and differential tt cross sections using the full data samples 
collected by the CDL and DO collaborations at the Tevatron is moving towards its completion. 


5 

























































Top Quark Pair Production Cross Section at the Tevatron 


Reinhild Yvonne Peters 



m(tt) [GeV] 


ra 

«} 

■o 

o 

4i4 

o 

'i? 

ra 


DC 


(b) D0 .i^ = 9.7fb-i 

- ■ Axi200L - Axi200A 

- ■ Axi2000L AXI400A 

— AxiSOOA 

AXI200R - Axi2000A 

— Axi2000R — Z’220 

■, 1 _L 



-T- 



. 

. L.. 


400 600 800 1000 1200 


m{ti) [GeV] 


Figure 3: The tt cross section measured as function of m„- and compared to various axigluon models 


These legacy measurements are important to test the SM, tune Monte Carlo simulations and search 
for hints of new physics. The different initial state and energy of the Tevatron compared to the 
LHC, make these complementary measurements. 
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